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ABSTRACT 


The metal matrix composites, have proved to be an important class of materials, 
with the potential to replace a number of conventional materials being used iu 
automotive, aerospace, defense, and leisure industries, where the demand for light weight 
and higher strength components is increasing. During centrifugal casting of Aluminium 
melts containing suspended ceramic particles, segregation of the particles occurs either 
towards the outer or towards the inner periphery of the casting depending on the relative 
density compared to that of the melt. A one dimensional heat transfer model coupled with 
equation for force balance on particles is developed to predict the temperature 
distribution in the casting and mold regions, solidification time of the casting and particle 
segregation pattern in the casting region. The model takes into consideration the 
propagation of solidification firont and movement of particles due to centrifugal 
acceleration which takes place either in the same or in opposite direction to that of the 
solidification front depending on the relative density difference between the particles and 
melt. In the force balance expression, repulsive force term is incorporated for the 
particles that are at the vicinity of the solid/liquid interface to calculate the particle 
segregation pattern in the casting region. The solution of the model equations has been 
obtained by pure implicit finite volume technique with modified variable time step 
(MVTS) approach. 

The effects of various process parameters such as, rotational speed of the mold, 
size of the reinforcing material, relative density difference between the particle and melt, 
initial pouring temperature of the liquid melt, mold pre-heating temperature, heat transfer 
coefficient between the casting/mold interface are studied. It is noted that for a given set 
of operating conditions, the thickness of the particle rich region in the composite 
decreases with increase in rotational speed, particle size, relative density difference 
between the particle and melt, initial pouring temperature and initial mold temperature. 
With decrease in the heat transfer coefficient between the casting/mold interface, the 
solidification time increases which, in turn, results in more intense segregation of solid 
particulates. Again, with increase m the initial volume fraction of the solid particulates, 
both the solidification time as well as the final thickness of the particulate rich region 
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increase. It is seen that, for AI-AI2O3 and Al-SiC systems, depending upon the variables 
used in the investigation, up to certain distances from the outer periphery of the casting, 
the volume fraction of the particulates in the solidified composite remains the same as 
that of the initial melt for all the three particle sizes considered here. But for Al-Gr 
system, this behavior is seen only in case of 1 pm particle size, for a particular set of 
parameters investigated here. 
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CHAPTER 1 


INTRODUCTION: 

1.1 INTRODUCTION: 


Composites can be considered materials consisting of two or more chemically distinct 
constituents, on a macro scale, having a distinct interface separating them [1]. Properties 
of composites are strongly influenced by the properties of their constituent materials, 
their distribution, and the interaction between them. These materials have many 
advantages over conventional materials because of then: superior performance in terms of 
strength, stiffiiess, toughness, high-temperature performance and other mechanical 
properties. Depending on the matrix materials, composites can be classified as poljmier 
matrix composites, metal matrix composites and ceramic matrix composites. 

Among this, metal matrix composites, where a strong ceramic reinforcement is 
incorporated into a metal matrix have proved to be an important class of materials, with 
the potential to replace a number of conventional materials being used in automotive, 
aerospace, defense, and leisure industries, where the demand for light weight and hi^er 
strength components is increasing. Metal matrix composites possess high levels of 
strength and stiffiiess, improved thermal stability, low coefficient of thermal expansion 
and improved wear and seizure resistance [2]. 

The functional gradient composite is a new type of composite that has become 
attractive because of its multi-functional properties, including the reinforced ceramic 
particle properties, matrix metal properties, and their combined properties [3]. These are 
designed to have heat resistant ceramics on their high temperature side and tou^ metal in 
their low temperature side, with a gradual compositional change. The main advantages of 
the functionally gradient composites are die gradual change of properties such as 
coefficient of thermal conductivity and electrical conductivity, thermal stress etc. 
resulting fi:om the gradual distribution of its structure and reinforced phases, which can 
avoid the destruction caused by the properties mismatch at high temperature. 
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This thesis deals primarily with the functionally graded metal matrix composites 
where the reinforcement is in the form of ceramic powders such as alumina, silicon 
carbide and graphite in aluminium matrix. 

1.2 DIFFERENT FABRICATION TECHNIQUES FOR THE PRODUCTION OF 
METAL MATRIX COMPOSITES: 

Ceramic reinforced metal matrix composites can be prepared by both solid state 
(including powder metallurgy and diffusion bonding) and liquid state (squeeze casting, 
pressure infiltration, stirring casting, centrifugal casting etc) [3] processes. However, 
solid-state processing techniques require more time, costly equipments and often these 
are cumbersome. Compared to solid-state processing methods, liquid state processing 
(like casting techniques) are characterized by several advantages, i.e., 

• Lower production cost 

• Time consumed is less 

• Simple in operation 

• Secondary processes are minimized 

• Intricate shapes can be produced 

Due to these advantages, liquid state processing techniques are becoming more 
popular. 

1.2.1 FABRICATION TECHNIQUES OF FUNCTIONALLY GRADIENT 
COMPOSITES: 

The fabrication of the functionally gradient composites has been accomplished 
via powder metallurgy (PM) techniques and deposition processes such as plasma 
spraying, physical vapor deposition (PVD), chemical vapor deposition (CVD), electro 
forming, and self - propagating high- temperature synthesis (SHS) [3]. Compositional 
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gradient is achieved by varying the chemistry of the powders and reactive gases with time 
in PM, PVD and CVD processes, respectively. But these processes are hard to be applied 
in practical use due to complex techniques and high cost. This has led to the development 
of an alternative technique. Centrifugal casting of metal matrix composites is an 
alternative technique to produce gradient composites very easily by controlling various 
process parameters according to the requirements. 

1.3 CENTRIFUGAL CASTING OF METAL MATRIX COMPOSITES: 

Centrifugal casting involves pouring liquid metal into a rapidly rotating mold, which 
may be moimted either vertically or horizontally, and continuing the rotation until 
solidification is complete. The main advantage of this technique is good mold filling 
combined with good micro stmctural control which usually gives excellent mechanical 
properties. During centrifugal casting of metal matrix composites, segregation of particles 
occurs due to centrifugal forces, either to the inner or to the outer part of the casting, 
depending on the relative densities of the particles and the melt, which results in particle 
reinforced functionally gradient composites [4-6]. The extent of segregation depends on 
various process parameters, such as geometry, pouring temperature, solidification time, 
density difference between the matrix and reinforced particles, rotational speed, etc. 

1.4 REQUIREMENTS OF MATHEMATICAL MODELLING OF 
CENTRIFUGAL CASTING OF METAL MATRIX COMPOSITES: 

During centrifugal casting of particles reinforced metal matrix composites, it is 
difficult to determine the temperature distribution and solidification time by experimental 
techniques. Because of this, accurate data on solidification time are not available. 
Estimation of solidification time and temperature distribution is a complex problem, 
because the moving ceramic particles disturb the temperature equilibrium. Also, the 
segregation of the particles alters the thermo-physical properties of composite melts 
during solidification. Therefore, it is desirable to estimate the solidification time and also 
influence of other processing parameters on the solidification time of metal matrix 
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composites through some alternate means. Mathematical modeling of the casting process 
based on heat and mass transfer analysis can be a useful alternate methodology. 

1.5 BRIEF REVIEW OF MATHEMATICAL MODELLING OF CENTRIFUGAL 
CASTING OF METAL MATRIX COMPOSITES: 

Kang and Rohatgi [7] have described the results of a heat transfer analysis of 
centrifugal casting of metal matrix composites by one-dimensional analysis considering 
the thermal mechanical properties due to particles moving as a function of temperature. 
These theoretical predictions have been compared with the results of Lajoye and Suery 
[8]. In their investigations, the positions of the dispersed particles at a given time have 
been analyzed as a first step. Then the temperature distributions in the mold and the 
solidifying metal have been analyzed at different time intervals and using these 
temperature distributions, solidification times for different centrifugal speeds, initial mold 
temperature, and pouring temperature of molten metal have been estimated. Siva Raju 
and Mehrotra [9] have presented a more realistic model where the volume fraction of the 
particles across the thickness of the casting varies with time. Siva Raju and Mehrotra 
have also considered variations in heat transfer coefficient and latent heat release. 

1.6 OBJECTIVE: 

The mam objectives of the present investigation are as follows: 

• To develop a model to describe the temperature distribution throughout the 
casting and mold regions, and particle segregation in the casting region. Finite 
Volume Technique has been used as the solution scheme unlike the finite 
difference method which was employed by Siva Raju and Mehrotra. 

• Incorporation of repulsive force term in the force balance expression to calculate 
the particle segregation pattern in the casting region. 

• To study the influence of different process parameters such as rotational speed of 
the mold, particle size, relative density difference between the particle and melt. 


4 



initial mold temperature, pouring temperature, heat transfer coefficient between 
the casting and mold on the solidification time of the casting and on the 
segregation pattern of the particles. 
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CHAPTER 2 

MODEL FORMULATION: 

Siva Raju and Mehrotra have presented a mathematical formulation based on one- 
dimensional heat transfer analysis incorporating variation in thermo-physical properties 
due to particle movement in the matrix. The present investigation is an attempt to further 
improve that analysis. At the outset, there are two principal differences between the 
present investigation and that of Siva Raju and Mehrotra. 

1. Both Kang and Rohatgi [7] and Siva Raju and Mehrotra [9] iu their models 
neglected the repulsive force term m their equations for force balance on particles. 
In the present investigations, the repulsive force which a particle experiences 
when it approaches a solid wall/solidification front is taken mto consideration and 
characterized by an appropriate expression, the details of which are given in a 
later section of this chapter. 

2. In the previous work of Siva Raju and Mehrotra, variation of volume fraction of 
particles with time is considered by ensuring conservation of total mass of 
particles at all times. In their work, the volume of particulates in between two 
consecutive nodal points always remain same. So, in that case the length of the 
segment varies with time and so are the nodes as they depend on constant volume 
fraction of particulates. But in the present formulation the nodes are fixed for all 
time and the variation of volume fraction of particles in any particular segment 
with time is calculated by taking into account the particle movement with time. 
One of the important features of the present model is the use of finite volume 
technique as the solution procedure unlike the finite difference scheme used in the 
previous work [10]. The advantage of finite volume method over finite difference 
one is the automatic conservation of all ph)^ical quantities unlike the 
approximation of so in the later case. The detailed description of the same is given 
in the chapter. 

A schematic representation of the centrifiigal casting of metal matrix 
composite, mathematically modeled in this investigation, is shown in Fig. 2.1. The 
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heat is withdrawn from the liquid region of the casting, which is at temperature of Tic , 
throu^ the solidified composite region at teiiq)erature Tsc, to graphite mold at 
temperature Tg and subsequently throu^ steel mold at temperature Tm to the 
surroundings. Heat is also radiated away from the iimer surface of the casting. The 
dispersed particles in the liquid are at temperature Tic. As the solidification proceeds by 
conduction of heat transfer through the composite to the graphite mold, the solid-liquid 
interface moves away from the graphite mold. 

2.1 MODEL FORMULATIONS TO KNOW THE TEMPERATURE 
DISTRIBUTION IN CASTING AND MOLD REGION AND SEGREGATION 
PATTERN OF THE PARTICLES: 

A mathematical model is formulated to predict the temperature distribution 
throughout the mold and casting regions and particle segregation within the Uquid region 
of the casting. 

The model is based on the following assumptions: 

2.1.1 MODEL ASSUMPTIONS: 

• The heat flow is purely one dimensional and perpendicular to the mold wall. 

• The mold is filled with liqmd metal cont aining solid particulates instantaneously, 
and that at the instant of mold fillin g the solid particulates are homogeneously 
distributed in the liquid metal matrix. 

• Thermal properties of solid particulates and metal matrix are temperature 
mvariant. 

• Heat transfer co-efficient between casting and graphite mold decreases as 
solidification proceeds due to increase in air gap because of contraction of 
casting. 

• Natural convection and movement of particles due to buoyancy are neglected. 

• The interface position between solid and liquid regions is calculated by assuming 
it to be planar. 
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• There is no thermal resistance between particles and liquid metal. 

• Reinforcement particles are assumed to be spherical in shape. 

• The motion of the particles is stopped by the liquidus front, i.e. the particles do 
not move in the mushy zone of the casting and they are not rejected by the 
solidification front. 

• For large volume fractions, interaction between particles occms thus reducing 
their velocity. It is assumed that the reduction in velocity can be characterized by 
an increase of apparent viscosity of the liquid. The apparent viscosity can be 
represented as [1 1] 

Vc= V [1+ 2.5 Vf (t) + 10.05 Vf^ (t)] (2.1) 

• Segregation of particles can only occur up to a maximum volume fraction equal to 
52 % [12]. 

2.1.2 HEAT CONDUCTION FORMULATION: 

This formulation is essentially the same as that of Kang and Rohatgi [7] and Siva 
Raju and Mehrotra [9]. 


2.1.2.1 Governing Equation: 

The heat transfer processes in various regions of casting and mold are governed 
by one dimensional unsteady state heat conduction equation written in cylindrical co- 
ordinates. 


Pi^i 


^=}.A.(rk 

dt r 5r ^ 


where. 


( 2 . 2 ) 
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z - Ic, sc, g, m for the liquid composite, solid composite, graphite and for steel mold 
regions, respectively. 


he , he , Csc , cic„ Psc. Pic are estimated by deterrnirdng the volume fractions of the 
solid particulates in the liquid and solid composite regions, which vary with time due to 
the movement of particles owing to the density difference between the liquid metal 
matrix and solid particulates. 

2.1. 2.2 Thermo-Physical Properties of Composites: 

Thermal conductivities, densities and specific heats of composites in either liquid 
or solid regions are determined using the rule of mixtures described as follows as a 
function of volume fraction of particulate, Vj(t) at any time ‘t’ [13]. 

^c=(i-v/f;;p„ + v/t;p^ (2.3) 

where. Pc refers to the property of composite under consideration, and Pp are the 
corresponding properties of the matrix and particulates , respectively. 

The volume fraction depends on the viscosity of the alloy, particle size, mold 
rotational speed and the density difference between the particle and the molten metal. 

2.1.2.3 Initial Conditions: 

During centrifugal casting, before pouring the molten metal (which is at 
temperature into the mold, the mold is preheated to a certain temperature (Tm) to 
avoid the thermal damage to the mold. Therefore the initial temperature distribution (at 
time /=0) in the casting and mold regions considered as 


Tic Tp 

(2.4) 

II 

II 

(2.5) 
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2.1.2.4 Boundary Conditions: 


The boundary conditions in different regions of the casting and mold are as follows: 


1 . At the inner surface of the casting, i.e. , at r = Ric 



dr 


= K(T,-T,) 


( 2 . 6 ) 


2. At the outer surface of the casting, i.e., at r = Roc 


-k.,^ = K(T^-Tgi) 
or 


(2.7) 


5. At the inner surface of the graphite mold, i.e., at r = Rig 


dT^ 


( 2 . 8 ) 


4. At the outer surface of the graphite mold, i.e., at r = Rog 


dT^ 

^ dr so miJ 


(2.9) 


5, At the inner surface of the steel mold, i.e., at r = Rim 


psn-i 


dr 


go 


■T„i) 


( 2 . 10 ) 


6 . At the outer surface of the steel mold, i.e., at r = Rom 
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( 2 . 11 ) 




7. At solid-liquid interface, i.e., at r = Rs(tj 

Tsc=T,c = Tf (2.12) 

8. The energy balance at solid-liquid interface, i.e., at r = Rs(t) is obtained by equating 
the rate of heat removed from the solid phase in the positive r direction to the sum of 
rate offbeat supplied to the interface from the liquid phase in the positive r direction 
and rate of heat liberated at the interface during solidification, i.e. 


, Sr, , dT ^^ds{t) 


SC ^ 

or 


dr 


dt 


(2.13) 


2.1. 2.5 Formulation of Heat Transfer Co-Efficient: 

The rate of solidification of the liquid composite is significantly dependent upon 
the air gap formed at the casting-graphite mold interface due to contraction of the casting 
as well as thermal expansion of the mold during solidification, and also to some extent on 
the air gap at graphite-steel mold interface due to imperfect contact. It is assumed that the 
heat transfer coefficient between casting and graphite mold due to air gap varies as [8]. 


h, 


= h, 




5(f) 


(2.14) 


where hi is initial heat transfer coefficient, hf is final heat transfer coefficient, s(t) is 
solidified thickness and u is total thickness of the casting. 
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2.1.3 FORCE BALANCE ON PARTICLES: 


Here it has been assumed that fluid is under laminar flow < I ). So, the force 
balance equation is as follows; 



liPp -Pi)a}^r~67rvRpV^ 



dt 


Under zero acceleration, the velocity of the spherical particle according to Stokes law is 


ARlKpco'^r 

ISv 


(2.15) 


where, Rp is radius of the particle, a> is angular velocity, r is particle position, v is 
viscosity of the melt and Kp = p^- i.e., density difference between the particle (pp) 

and liquid metal (pi) .The Reynolds number is given by 


R. = 


PiVAR, 


(2.16) 


During vertical centrifugal casting, a particle which is suspended in the liquid is 
subjected to a vertical acceleration due to gravity g and a centrifugal acceleration 
y = co^r. Generally y is much greater than g which allows the vertical displacement of 
the particle to be ignored. Therefore the different forces on a particle are 

• Centrifugal force due to tiie rotation 

• Viscous force due to drag effect 

• Repulsive force due to movement of solid-liquid interface 

A schematic representation of the various forces acting on the particle is given in Fig. 2.2. 
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SOLID 


ds{t) 

V = — 

dt 

LIQUID 




Figure 2.2: Schematic diagram of various forces acting on a moving 
particle in the liquid melt 



The force balance equation on the particle due to centrifugal, viscous, and repulsive 
forces can be expressed as [14] 

F^-F,-FR=Fnet (2.17) 

where, F„et is the net force on the particle. Fa, is the force due to centrifugal acceleration, 
Fv is viscous force, and Fr is the repulsive force. It is to be noted that the repulsive force 
is significant only on those particles which are in the vicinity of a solid wall or the solid- 
liquid interface. Hence, the force balance for the particles which are far away from the 
solid-liquid interface can be written as follows. 

F^ -Fy=F„e, (2.18) 

iPp - PiWr - SnvRp ^ = ^^IPp ^ (2-19) 

Solution of Eq. (2. 19) for a particle moving at a constant velocity, gives its position at 
any instant of time t as 


r(t) = exp 


4o}\pp-p,)Rlt 

18v 


( 2 . 20 ) 


where ro(t) is the position of the particle at time t=0 
2.1.3.1 Repulsive Force : 

2.I.3.I.I. Importance of the use of repulsive force term: 

In the present work of centrifugal casting of metal matrix composites, the final 
distribution of the second phase in the solidified material is important for obtaining the 
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requisite properties, hence an understanding of the phenomena which contribute to the 
redistribution of the particles during solidification is very important in order to achieve 
the desired results. 

When the solidification front of an alloy approaches a particle in the melt, 
interactions take place between them due to three main reasons: the molecular surface 
forces between the particles and the front, distortion of the temperature field because of 
the different thermal properties of the particle and distortion of the solute concentration 
field due to blocking of solute diffusion from the front by the particle. If the molecular 
surface forces are repulsive, the front can push the particles along with it unless the 
repulsive force is overcome by the viscous drag on the particle [15]. Several theories 
have been proposed which are based on the concept of molecular surface forces repelling 
the particle from the growth front and thus facilitating the entry of liquid into the gap 
between the growing solid and the particle [16-18]. The magnitude of the forces depend 
on the dimensions of the gap between the particle and the front and hence on the 
morphology of the front. When a liquid containing a dispersion of particles is solidified, 
the solidification front can exert a repulsive force on the particles and push them along 
with it, thereby increasing the particle concentration in the last solidifying liquid. The 
pushing of the particle is considered as a steady state phenomenon in which the particle 
and the solidification front move with the same velocity and the gap width is maintained 
constant. 


2.I.3.I.2. Proposed Theories for the Expression of Repulsive Force: 

Several models explaining the concept of repulsive force or wall effect on the 
particle or individual have been developed when the laters enter into the influence zone 
of the wall or solid surface. 

The model proposed by Sasikumar and Kumar [19] demonstrates a system 
consisting of silicon carbide particles dispersed in Aluminium melt, solidified vertically 
from the bottom as well as in cylindrical molds. They have explained the effect of 
repulsive force on the particle exerted by the soUdification front. It has also been known 
that the phenomenon of particle pushing occurs only below a certain velocity of the 
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advancing solidification front known as the critical velocity. If the velocity of the growth 
front is greater than the critical velocity, the particle is engulfed by the front. The critical 
velocity is a function of the particle size, the temperature gradient at the solidification 
front, and a number of material properties like the ratio of thermal conductivities of the 
particle and the melt, the co-efficient of molecular surface forces between the particle and 
the solidification front, etc. In this method for calculating the critical velocity, the first 
step is the calculation of the change m shape of the solidification front as it approaches 
the particle. From this, the dimensions of the gap between the particle and the front are 
calculated. The repulsive force sxsrtsd. oil the particle depend on the gap and these are 
calculated and balanced to obtain the velocity of the particle. 






where. 



( 2 . 21 ) 


( 2 . 22 ) 


rp = Radius of the particle 

h(r) = The gap between the particle and the front 

B = The co-efficient for molecular interaction between the front and the particle 
T] = Viscosity of the melt 

The maximum velocity under which steady state is possible is identified as the critical 
velocity. Numerically this is determined by calculating the Vp at different distances, L of 
die solidification front from the particle. As L decreases, Vp mcreases and then goes 
throu^ a maximum. It is concluded from this that if the solidification front were to have 
a velocity, V, equal to Vp(L) , steady state pushing of the particle occurs when the front is 
at a distance L from the particle. If the V were to exceed the maximum possible value for 
Vp then no steady state pushing is possible; the gap between the particle and the front 
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continues to decrease and the particle is engulfed by the front. The maximum value of Vp 
is therefore identified as the critical velocity. 

Sannomiya and Matuda [20] has proposed a model which describes the behavior 
of fish in a water tank where the wall exerts the repulsive force on the individual when it 
approaches towards it. The detailed description of this model is given in the later section 
of this chapter. 

In this thesis, the expression of the repulsive force that is exerted by the 
solidification front on the particle is taken from the mathematical model that described by 
Sannomiya and Mafiida [20]. 

2.I.3.I.3. Reasons in support of the use of the same expression of the repulsive 
force in the case of a particle: 

As fish is a living being, its movement is influenced by several factors such as, its 
velocity relative to the ambient water, the extent of swimming forward continuously 
according to its swimming ability, interaction among the individuals, the character of 
uniforming the swimming speed and the direction in the school, repiilsion and attraction 
from the walls when the individual approaches towards it closely. The forces which 
describe all these activities are damping force, propulsive force, interactive force, 
schooling force, repulsive and attractive force from the wall respectively [20]. 

As fish is a living species, when it approaches the wall, it avoids to strike against 
the same which is a psychological phenomenon happening to all living beings. But in the 
case of solidification of a liquid containing a dispersion of second phase particles, when 
the later reaches the solidification front, they do not re-bounce (due to the involvement of 
no psychological phenomenon and the assumption that there is no rejection of particles 
by the solidification front ). The velocity of the particles slow down in the influence zone 
of the wall due to wall effect and finally they get absorbed by the solidification front. As 
repulsive force is independent of the tank shape and size [21], so the parameters 
estimated by Sannomiya and Matuda is used in the present thesis. 
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2.1.3.1.4 Expression for repulsive force : 


For the particles which are in the vicinity of a solid surface or 
interface, the expression of the repulsive force is [20] 

/ = 1 


/ 


wil"^ 



10 


forv„ > Oanddj, < d* 

otherwise 


where is a constant co-efficient . 

L is the number of wall sides. 

is the width of the influence zone of the wall (« 10'® cm) [22] 
dll is the distance from the individual z to the wall/, 
the unit vector e/ is normal to the wall. 

For the individual i, vu is the velocity component normal to the wall /. 
So, particle position at any given instant of time is given by 


r(0 - ro exp 


4/3nGi^{pp-pl)Rpt 

K+6kvRp 


where, K = k ^ 


d* 


the solid-liquid 


( 2.23) 


(2.24) 


(2.25) 


( 2.26) 
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2.1.3.1.5 Estimation of by Sannomiya and Matuda: 


In the work of Sannomiya and Matuda, an individual fish is assumed to be a 
particle. The depth of water in the tank where experiment has been carried out is shallow 
and the motion of fish is restricted within a two-dimensional space. Since only a small 
school is considered, the motion of each fish is described individually by a state equation. 
The experiment used to start with releasing the fish into a specific position of the tank 
and the motion of fish is restricted within a two-dimensional space. As a general rule, the 
fish school swims along the wall. The image of the fish behavior used to be recorded 
during ten minutes by a video tape recorder. The image data was sampled at intervals of 
O.Ssecs to obtain the time series data. Since only a small school has been considered , the 
motion of each fish is described individually by a state equation. The parameters included 
in the equation are identical irrespective of the individual, under the assumption that fish 
in a school are homogeneous. When an individual approaches a waU closely, it avoids 
striking against the wall. Repulsive force expresses such a behavior. Sannomiya and 
Matuda have estimated the value of kj" by using the observation data of this water tank 
experiment and by applying the least squares algorithm. The validity is investigated by 
testing the whiteness of the residual of the state equation and also by comparing the 
simulation result with those of the experimental ones. 

2.1.4 DETERMINATION OF VOLUME FRACTION : 

In order to determine the volume fraction variation in liquid metal matrix, the 
casting thickness is divided into n equal segments. Nodal pomts are considered as the 
extreme positions of each segment. The volume fraction of particulates in a particular 
segment is defined as the ratio of volume of particles to the total volume of that segment. 
This is given as: 


Vfs = 


Vs 

Vs + Vi 


1 


1 + 


psmi 


pirns 


(2.27) 
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where is the volume fraction of the particulates, Fs is the volume of reinforced 
particles, and F/ is the volume of the metal matrix in each segment. Since the initial 
volume fraction of the particulates is known in each segment, the volumes and masses of 
the solid particles and matrix melt can be obtained using the following relationships. 

Fs=FfsF 

Fi = F- Fs (2.28) 

m,== V,p, (2.29) 

where F is total volume of each segment, m, is mass of solid in each segment, and mi is 
total mass of liquid in each segment. 

For simplicity, initially the particle positions are considered at the nodal 
points at time t~0. The particles at different nodal points have different velocities, 
because the motion of each particle in the liquid melt is dependent on its position. The 
total number of particulates in the casting region always remains the same. The new 
particle positions for time t+A are obtained by using either Eq. (2.20) or (2.25) taking 
into account the appropriate conditions. The new particles may fall either on any node or 
in the interior of any volmne segment. From this, the total number of particulates in any 
volume segment can be found out and thus, the new volume of the particulates in each 
segment can be calculated. The volume of each segment remains fixed for all time. So 
from this, the new volume fraction of the particulates per unit length can be obtained for 
time t+At. To calculate the particle position for the next time step, initial positions of the 
particulates are updated with the new particle positions tiiat have already obtained and 
further calculations are carried out in the similar manner. 
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CHAPTERS 


SOLUTION PROCEDURE; 

The model equations have been solved numerically by using simple implicit finite 
volume technique. For this the “ r - r ” domain is subdivided into small intervals of 
constant Ar in space and variable At in time as shown in Fig. 3.1. The variable time step 
approach is used to solve the problem. This approach requires that at each time level 
the time step At n is so chosen that the interface moves exactly a distance Ar during the 
time interval At, hence always stays on the node. Therefore, the problem is mainly 
concerned with the determination of the time step At= t„+; - t„ such that in the time 
interval from t„ to t„+;, the interface moves from the position nAr to the next position 
(n+l)Ar[23]. 


t 



0 lAr 2Ar 3Ar nAr (n+l)Ar 


Figure 3.1: Subdivision of“r-t” domain using constant Ar, variable At 
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3.1 ADVANTAGES OF FVM OVER FDM: 


The usual procedure for deriving finite-difference equations, consists of 
approximating the derivatives in the differential equation via a truncated Taylor series. 
The method includes the assmnption that the variation of the dependent variable is 
somewhat like a polynomial in the independent variable, so that the higher derivatives are 
unimportant. This assumption, however, leads to an undesirable formulation when, 
exponential variations are encountered. The Taylor-series formulation is relatively 
straight-forward but allows less flexibility and provides little insist into the physical 
meanings of the terms [24]. 

But, the basic idea of finite volume formulation lends itself to direct physical 
interpretation. Here, the calculation domain is divided into a number of non-overlapping 
control volumes such that there is one control volume surrounding each grid point. The 
differential equation is integrated over each control volume. Piecewise profiles 
expressing the variations of the dependent variable between the grid points are used to 
evaluate the required integrals. The result is the discretization equation containing the 
values of dependent variables for a group of grid points. The discretization equation 
obtained in this manner expresses the conservation principle for the dependent variable 
for tbe finite control volume, just as the dififerential equation expresses it for an 
infinitesimal control volume. The most attractive feature of the control-volume 
formulation is that the resulting solution would imply that the integral conservation of 
quantities such as mass, momentum, and energy is exactly satisfied over any group of 
control volumes and, of course, over the whole calculation domain. This characterstic 
exists for any number of grid points, not just in a limiting sense when tbie number of grid 
points becomes large. Thus, even the coarse-grid solution exhibits exact integral balances 
[24,25]. 
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3.2 VALIDATION OF FINITE VOLUME TECHNIQUE: 


The validation of this technique has been carried out by comparing the numerical 
solutions with those of the analytical ones in determining the temperature profiles for 
various cases. Figures 3.2 and 3.3 show the comparison between the analytical and 
numerical techniques for unsteady-state heat conduction in a slab and cylinder of finite 
thicknesses, where the initial temperatures are given and surface temperatures are fixed 
[26-29] and Fig. 3.4 shows the validation of this technique for the steady state 
temperature profiles in the case of a finite hollow cylinder where temperature is fixed at 
one boundary and the other boundary is a convective one. The analytical solution for Fig. 
3.4 is as follows: 


hr, . 


hr. 

frA 

1-1- —Lin 


T = T, +— Lr,hi 

2 

k 

U j 

* 



(3.1) 


where, 

T 2 = Temperature of the outer surface of the cylinder, i.e., at r = 

Tf= Temperature of the ambient adjacent to the inner surface of the cylinder 
h = Heat transfer coefficient at the inner surface of casting 
k = Thermal conductivity of the cylinder 
r = Distance firom the axis of the cylinder 

3.3 FINITE VOLUME APPROXIMATION; 


The differential equation and the boundary conditions for both the mold and 
casting regions can be discretized by using implicit method. The thickness of each mold 
region is subdivided into n equal grids. The number of grids in the solid and liquid 
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Figure 3.2: Temperature profiles for unsteady-state heat conduction in a slab of finite 
Thickness 
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Figure 3.3: Temperature profiles for unsteady-state heat conduction in a cylinder of finite 
Thickness 
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Figure 3.4: Temperature profiles for steady-state heat conduction in a cylinder of finite 
thickness 
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composite regions vary with time. It means that for the solid composite region the 
number of grids («i) goes on increasing as the solidification proceeds, and for the liquid 
composite region, the number of grids (n-ni) decreases, but their sum always remains 
same as n. 

Equation (2.2) is solved with boundary conditions Eqs. (2.6) through (2.12) using 
pure implicit finite volume technique. The resultant discretized equations are arranged in 
a tri-diagonol matrix form and the solution of these equations is obtained by using 
Thomas Algorithm (TDMA) v/hich gives temperature distribution in both casting and 
mold regions for a particular time step Al„_ The detailed description is presented in 
section 3.5 of the present chapter. 

3.4 DETERMINATION OF TIME STEPS: 

During solidification, the interface moves from the position n A r to the position 

(n+i) Ar within a time interval of At ,= t/+; - 1,-, which is not known in a priori. 

Using modified variable time step (MVTS) method, the actual value of At /is obtained by 
iteration as follows [30]. 

1 . An initial guess value for the time step At is chosen as At i = At P . 

2. Using At file nodal temperatures are obtained by solving the finite volume equations 

for the mold and casting regions. 

3. Then the time step At is obtained by using the estimated temperatures in step (2) in 
Eq. (3.13). 

4. Using At , the steps (2) and (3) are repeated to obtain At . 

5. Iteration is carried out until the difference between the two consecutive time steps 
satisfies a specified convergence criterion. 

6. To obtain A t,+/ , At|-^|* = At is assumed and the steps from (2) to (5) are repeated. 

7. The steps from (2) to (6) are repeated till the solid-liquid interface reaches the inner 
surface of the casting. 
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3.5 DISCRETIZATION OF THE GDE AND BOUNDARY CONDITIONS 
ACCORDING TO THE CONTROL VOLUME TECHNIQUE: 



Figure 3.5: Schematic diagram showing different nodal points of the casting and mold 
Region 
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“Half’ control volume 



Figure 3.6: Control volumes for the internal and boundary points 


FOR ALL INTERNAL POINTS: 



r 


Figure 3.7: Grid- point cluster for the one-dimensional problem 
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The differential equation (2.2) is integrated over the full control volume as follows: 


J ]rfx:~dtdr= J \^(rkf)drdt 


Then the discretized form of Eq. (3.2) becomes: 


( rk\( rk 




+ tH- + 7^^ n = h + T- 


FOR INNER SURFACE OF THE CASTING ( NODE ‘0’): 


Integrating the differential equation (2.2) over half control volume: 


j Jrpc—dtdr= J j—{rkf)drdt 


and substituting the boundary condition [Eq. (2.6)], the discretized equation becomes: 


29 / [dr J 




where. 


^ _ 'Ik-J^Q 
V^i+^o 


= thermal conductivity of node ‘ 1 ’ 
kQ= thermal conductivity of node ‘0’ 
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FOR NODE ‘a’: 


Integrating the governing equation (2.2) over the full control volume and using the 
boundary condition [Eq. (2.7)], the discretized equation becomes: 


(r+r) 

pcdr^^^—^ 1 + 


2dt 


2k 




sc J 


iPr). 




'{T = 


^rXhP 

e e \ 


2k 
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sc J 


^ r k ^ 

{Sri 


T 

^a-l 



(3.7) 

FOR NODE ‘b’: 


Integrating the governing equation (2.2) over the full control volume and using the 
boundary condition [Eq. (2.8)], the discretized equation becomes: 


pcdr 


{re + rj 
2dt 
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(3.8) 


FOR NODE ‘c’: 

Integrating the governing equation (2.2) over the full control volume and using the 
boundary condition [Eq. (2.9)], the discretized equation becomes: 


r 


.{r.+rj 


pcdr- 
\ 2dt 
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(dr) 


yj 


^ r,kX ^ 
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2dt 


(3.9) 
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FOR NODE ‘d’: 


Integrating the governing Eq. (2.2) over the full control volume and using the boundary 
condition [Eq. (2.10)], the discretization equation becomes: 


If (r +r Y] 

^3.. V e ^ 'w/ 

1 


1 


T — 

(rk ^ 

e e 

T 4- 

f '] 

T 4- 


Aut j 
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■^d+\ ^ 

1 2A:, j 


2dt J 


FOR OUTER SURFACE OF THE STEEL MOLD ( NODE ‘e’ ): 

Integrating the governing equation (2.2) over the half control volume: 


e t+At 


I Jrpc^dtdr^ j 


r+A/ e 


dt 


idr 


(3.11) 


and using the boundary condition [Eq. (2.1 1)] , the discretized equation becomes: 


2dt 






+ 


Pe^p-e -P fii'^l 
Idt 


.2\ 


(3.12) 

At Solid-Liquid Interface: 

Equation (2.12) represents the condition of continuity of temperature at the interface, i.e.: 


Tsc =T,c = Tf 


( 3 . 13 ) 



The time required to move the solidification front from one node to the other can be 
derived from the Eq. (2.13) as follows: 


p.hM) 

(dT,\ , (dTA 
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sc I 
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1 J 


(3.14) 


Temperature profile for all the internal points, when the solidification front moves from 
one node to the other, is shown in Fig. 3.7. So, the time required to move the 
solidification front from one node to the other for all internal points is given as: 


At =• 






(M 


-k, 


sq 




t^-tA 

, (a^), 




(3.15) 


When the solidification front moves to frie last node, i.e., node adjacent to the ambient, 
then time required is: 


At =- 


Psc^e^it) 


^2 (^0 ) ^ SC ( 


r,. -Tf 




C J 


(3.16) 


and the temperature profile is shown in Fig. 3.8. 
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Figure 3.7: Schematic diagram of the temperature profile for all internal points 
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Figure 3.8: Schematic diagram of the temperature profile when the 

solidification front moves to the node adjacent to the ambient 
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Equations (3.3), (3.5), (3.7), (3.8), (3.9), (3.10), (3.12), (3.13) are arranged in the tri- 
diagonal matrix form and the solution of these equations is obtained by using the Thomas 
Algorithm (TDMA). This results in temperature distribution in both the casting and mold 
regions for a particular time step At . 

3.6 SEGREGATION OF PARTICLES: 

During centrifugal casting, segregation of particles takes place in the liquid 
compcsite due to the movement of particles resulting from difference in densities of tb.e 
particles and the melt, and also due to centrifugal acceleration. By solving Eqs. (2.27) and 
(2.28) using tbe appropriate equations from Eqs. (2.20) and (2.25) through Eq. (2.1) the 
thickness of the particle-rich region is estimated for various rotational speeds of the mold 
at various times. The particle movement is neglected in the mushy zone of the casting, 
which is characterized by the temperature region (Ts<T^Tl), because on this zone the 
viscosity, Vc, sharply increases. 
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CHAPTER 4 


RESULTS AND DISCUSSIONS: 

In this thesis, the model equations are solved to evaluate the followings for various 
operating conditions. 

i. Particle segregation pattern 

ii. Solidification time 

iii. Temperature distribution in the casting and mold region 
The parameters examined are: 

i. Rotational speed of the mold 

ii. Size of the reinforcing material 

iii. Relative density difference between the matrix and reinforcing material 

iv. Heat transfer co-efficient between the casting/mold interface 

V. Initial pouring temperature of the liquid melt 

vi. Initial mold temperature 

vii. Initial volume fraction of the particulates in the melt 

Thermo-physical properties of the matrix metal, reinforcement materials, graphite and 
steel molds used in this simulation are given in Table 4.1 [31,32]. Various design and 
operating parameters, like geometric constants for the casting and the mold, the heat 
transfer coefficients at different regions of the casting and the mold, and initial 
temperatures of the mold and metal used in the simulation are given in Table 4.2 [7,33]. 
The cooling conditions at the iimer surface of the casting and at the outer surface of the 
steel mold are defined in terms of the heat transfer co-efficient h 2 , hs respectively and the 
heat transfer due to the air gap at the metal-mold interface is characterized by hj. Since 
the actual values of hi are not known, computations have been carried out for a wide 
range of values from 1000 to 10000 W/m^/K for this co-efficient. The results of this 
simulation are presented in the following sections. 
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Table 4.1 : Thermo-physical properties of matrix metal, reinforcement particles, and 
mold materials [31,32] 


Thermo-physical 

Properties 

Aluminium 

(A356) 

AI 2 O 3 

SiC 

Graphite 

Carbon 

Steel Mold 

Graphite 

Mold 

k{Wm-^K‘) 

159 

24 

24 

38 

57.8 

38 

p(kgm'^) 

2685 

4000 

3200 

1900 

7800 

1900 

C(Jkg-‘K-’) 

963 

600 

690 

710 

481 

710 

TsCQ 

555 

- 

- 

- 

- 

- 

TlCQ 

615 

- 

- 

- 

- 

- 

TfCC) 

555 

- 

- 

- 

- 

- 

H(KJkg-^) 

389 

- 

- 

- 

- 

- 

v(Ns m'^) 

0.002 

- 

- 

- 

- 

- 
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Table 4.2: Design & operating parameters used in simulation [7,33] 


Outer diameter of steel mold, Dom (mm) 

215.0 

Outer diameter of graphite mold. Dog ( mm) 

150.0 

Outer diameter of casting. Doc ( mm) 

100.0 

Inner diameter of casting, Ac ( mm) 

80.8 

Initial heat transfer co-efficient due to air gap at the 
metal-mold interface, hi {W ) 

1000 to 10000 

Heat transfer co-efficient at the inner surface of 
casting, h 2 (JV m '^ ) 

8.4 

Heat transfer co-efficient at the outer surface of steel 
mold, h 3 {Wm-^K ^ ) 

8.4 

Heat transfer co-efficient between the graphite-steel 
mold interface, h 4 {W m'^K ^ ) 

10000 

Initial Pouring temperature of liquid metal, Tp (°C) 

650,730,800 

Initial mold temperature, Tm (°C) 

200,250,300 

Speed of rotation, N (rpm) 

1000,1500,2000 

Particle size, Rp (pm) 

1,2,3 

Ratio of initial to final heat transfer co-efficient, — 

10 





4.1 CALCULATION OF PARTICLE SEGREGATION PATTERN: 


The results of this simulation are presented in the form of 

i. Volume fraction of the particles 

ii. Thickness of the particle rich region for the entire region of the casting 
4.1.1 Effect of Rotational Speed of the Mold: 

Figures 4.1 and 4.2 show the effect of rotational speed of the mold on the 
segregation pattern of the particles, i.e., on the thickness of the particle rich region. The 
various rotational speeds used in these simulations are 1000,1500 and 2000 rpm 
respectively. Figure 4.1 is for the condition where the reinforcement particles are heavier 
than the liquid melt and Fig. 4.2 is for where the matrix metal is heavier than the 
reinforcing particles. From these two figures it is evident that, thickness of the particle 
segregated region decreases with increase in the rotational speed of the mold, i.e., more 
denser segregation pattern is observed with higher rotational speed. This is because, 
increase in rotational speed of the mold increases the centrifugal acceleration, which 
subsequently increases the centrifugal force on the particles. Depending upon the relative 
density difference between the particles and the matrix melt, and whether it is positive or 
negative, the particles move with a greater velocity towards the outer or inner periphery 
of the casting and get segregated there. It is observed in Fig. 4.1 that, there is a horizontal 
portion near the outer periphery of the casting which signifies that the final volume 
fraction of the particulates in the solidified composite is equal to that of the mitial volinne 
fraction in the melt, before reaching the maximum in the segregation pattern and this 
horizontal portion decreases with increase in the rotational speed of the mold. The reason 
is, rate of solidification is very fast at the beginning, which is more than the particle 
velocity due to which the particles are not able to move towards the casting/mold 
interface leading to the volume fraction of the particulates in the sohdified composite 
same as that of the initial volume fractions up to certain distances from the outer 
periphery. Gradually as the solidification front moves to the interior of the casting, time 
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complete solidification for Pp> pi 
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required to solidify a volume segment having the same dimension increases, as a result of 
which particles get more time to move towards the outer periphery of the casting leading 
to a segregated region next to the horizontal pattern. But with increase in the rotational 
speed of the mold, the centrifugal force on the particles increases and after certain time 
the velocity of the particles due to the centrifugal force becomes more than the velocity 
of the solidification front leading to less horizontal portion in the case of higher rotational 
speed. It is observed that, the thickness of the horizontal portion in case of Al- AI 2 O 3 
decreases from 9.73x10'^ to 5.21x10'^ m when the rotational speed increases from 1000 
to 2000 rpm. For Figure 4.2, particles get sufficient time to move to the inner surface of 
the casting, which leads to a maximum segregation at the inner periphery of the casting. 
The thickness of the AI 2 O 3 rich region varies from 3.085x10’^ to 2.337x10'^ m, while that 
of graphite rich region varies from 2.301x10'^ to 2.128x10'^ m as the rotational speed 
increases from 1000 to 2000 rpm after complete sohdification, as shown in Figs. 4.1 and 
4.2. 


4. 1 .2 Effect of Particle Size: 

Figures 4.3 and 4.4 show the effect of particle size on the thickness of the particle 
rich region, for three different particle sizes, namely 1, 2 and 3 /jm. Figure 4.3 is the case 
where the density of the particle is more than that of the matrix material and Fig. 4.4 is 
the case where the density of the particle is less than that of the liquid melt. From these 
two figures it is evident that, the thickness of the particle rich region decreases with 
increase m particle size leading to more intense particle segregation. With increase in 
particle size, both the centrifugal and viscous forces increase, but the increase in 
centrifugal force is more than that of the viscous force. Hence, larger particles move 
rapidly towards the outer or inner periphery of the casting depending upon the relative 
density difference between the particle and the liquid melt and form more intense 
segregation. In Fig. 4.3, the horizontal nature of the segregation pattern decreases with 
increase in the particle size. This is because, larger the particle size, larger will be its 
momentum. As a result, particles with greater momentum move faster than the rate of 
solidification front so penetrate into the melt to a farther distance than the particles of 
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Figure 4.3; Effect of particle size on particle segregation pattern after complete 
solidification for Pp > pi 



Figure 4.4: Effect of particle size on particle segregation pattern after complete 
solidification for Pp<pi 
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smaller size. It is observed that, the thickness of the horizontal portion in case of Al- 
AI 2 OS decreases from 1.933x10 to 6.85x10 m when the particle size increases from 1 
to 3 /jm. In Al-Gr system, for 1 ^ particle size the volume fraction of the particles in the 
solidified composite, in the volume segment just adjacent to the casting/mold interface 
remains the same as that of the initial volume fraction in the melt unlike to that of 2 and 3 
jiim particle size. The distance of this horizontal region from the casting/mold interface is 
3.4273x10’® m. The reason is, with increase in particle size, centrifugal force on the 
particle increases and larger particles move with hi^er velocity towards the inner 
periphery and do not get trapped by the solidifrcation front near the outer periphery of the 
casting. But in case of 1 fun particle size, the velocity of the particles in volume segment 
adjacent to the casting/mold interface is lesser than the velocity of the solidification front 
which results in trapping of the particles within the solidifying layer. As a result, in this 
volume segment the volume fraction of the particulates in composite is the same as that 
of the initial melt. The thickness of the AI 2 O 3 rich region varies from 4.449x10'^ to 
2.639x10'^ m, while that of graphite rich region varies from 2.519x10’^ to 2.259x10’^ m 
as the particle size increases from 1 to 3 after complete solidification, as shown in 
Figs. 4.3 and 4.4. 

4.1.3 Effect of Relative Density Difference Between the Particle and melt: 

Figure 4.5 shows the effect of relative density difference between the particles 
and the melt on the thickness of the particle rich region for the case where the density of 
the particle is more than that of the melt. It is seen that the thickness of the particle rich 
region decreases with increase in the relative density difference between the particle and 
the melt. The reason that with increase in the relative density difference between the 
particle and the melt, the net centrifugal force on the particle increases due to which the 
particles travel to a farther distance leading to more intense particle segregation. As the 
particulate density decreases from 4000 kg/m^ (AI 2 O 3 ) to 3200 kg/m^ (SiC), the thickness 
of the particle segregated region increases from 2.989x10 ^ to 4.003x10 ^ m. 


44 



4.1.4 Effect of Initial Pouring Temperature: 


Figures 4.6 to 4.8 show the effect of initial pouring temperature of the liquid melt 
on the particle segregation pattern for three different particle sizes, namely 1, 2 and 3 ism. 
The three different values of 7}, considered are, 650, 730 and 800 °C. All other parameters 
are kept fixed in these simulations. Figure 4.9 is a plot of variation of thickness of the 
particle rich region with initial pouring temperature for different particle sizes. It is 
observed that with increase in initial pouring temperature, the thickness of the particle 
rich region decreases for all the three particle sizes -that are considered. Higher initial 
pouring temperature essentially means that the larger amount of heat is to be removed 
from the melt before solidification begins. Hence, the metal matrix composites takes 
longer time to solidify and during this extra time the reinforcement particles are able to 
segregate more forming a denser segregation pattern. 



Figure 4.5: Effect of relative density difference between the reinforcement material and 
the melt on particle segregation pattern after complete solidification 
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Figure 4.7: Effect of pouring temperature on the particle segregation pattern for 2pm 
particle size 
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Figure 4.8: Effect of pouring temperature on the particle segregation pattern for 3/rm 
particle size 



Figure 4.9: Influence of initial pouring temperature (Tp) on the variation of thickness of 
particle rich region of Al-AhO 3 composite 
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4.1.5 Effect of Initial Mold Temperature: 


Figures 4.10 to 4.12 show the effect of initial mold temperature on the particle 
segregation pattern for three different values oiTu, namely 200, 250 and 300 °C and for 
three different particle sizes, viz 1 , 2 and 3 fjm. All the other parameters are kept fixed. 
Figure 4.13 is a plot of thickness of the particle rich region verses the initial mold 
temperature for the three different particle sizes. It is observed that with increase in initial 
mold temperature, the thickness of the particle rich region decreases as in the previous 
case. Higher initial mold temperature essentially means, reduced thermal gradient 
between the melt and the mold resulting in reduced heat transfer between them. Hence, 
the metal matrix composites take longer time to solidify and during this extra time the 
reinforcement particles are able to segregate more. 

4.1.6 Effect of Heat Transfer Co-efficient at the Metal/Mold Interface: 

Figures 4.14 to 4.16 show the effect of heat transfer coefficient at the metal/mold 
interface on the segregation pattern of the particles. It is seen that the thickness of the 
particle rich region increases with increase in heat transfer coefficient between the metal 
and the mold. This is because the rate of heat transfer is directly proportional to the heat 
transfer coefficient at the metal/mold interface. Higher the heat transfer coefficient at this 
interface, higher is the rate of heat transfer. Thus reduces the solidification time of the 
metal matrix and leads to an increase in particle rich region because the particles do not 
get enough time for more intense segregation and therefore remain distributed in the 
larger region of the composite. 

4.1.7 Effect of Particulate Volume Fraction: 

Figure 4.17 shows the effect of initial volume firaction of the particulates on the 
final distribution of the particles in the solidified composite and also on the thickness of 
the particulate rich region. It is evident from the figure that the thickness of the particle 
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Figure 4.10: Effect of mold temperature on the particle segregation pattern for l/jm 
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Figure 4.1 1 : Effect of mold temperature on the particle segregation pattern for 2/jm 
particle size 
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Figure 4.12; Effect of mold temperature on the particle segregation pattern for 3/jm 
particle size 


cji 


d 

o 

*5jd 

u 

u 

d 

a 


IZ> 

1/3 

0^ 

d 

o 

2 

H 


9.6 1 



Material : AI-AI 2 O 3 

- 

N : 1000 rpm 


hi : 1000 W/m*/K 

7.2 - 


4.8 - 

n- 


■° □ 


* ^ 

2.4 - 

0 ^ ^ ^ 

■ — -0 


— R = 1 urn 

- 

— R*" = 2 pm 

0.0 - 

— 0 — Rp = 3 pm 

1 1 1 1 1 1 1 1 1 1 1 1 1 


180 200 220 240 260 280 300 320 


Initial mold temperature (**C) 


4.13: Influence of initial mold temperature (Tm) on the variation of thickness of particle 
rich region of AUAhOs composite 
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Figure 4.15; Effect of heat transfer coefficient on the particle segregation pattern for 2 pim 
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Figure 4.16: Effect of heat transfer coefficient on the particle segregation pattern for 3 /ctj 
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Figure 4.17: Effect of initial volume fraction of particulates on particle segregation 
pattern 
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rich regions increases with increase in the initial volume fraction of the particulates. This 
is because, the increase in the initial volume fraction of the particles in the melt increases 
its viscosity, thereby reducing the motion of the particles which results in larger particle 
rich region due to reduced segregation. 

4.1.8 Effect of Solidification Time: 

Figures 4.18 to 4.23 show the variation of volume fraction of the particulates with 
time for three different particle sizes, namely, 1, 2 and 2>fjm. The heat transfer coefficient 
at the metal/mold interface is taken as 2000 W/m^/K in Figs. 4.18 to 4.20 and 5000 
W/m^/K for Figs. 4.21 to 4.23. All other parameters are kept fixed. From these figures it is 
evident that with increase in solidification time particles get more time to move in the 
liquid melt leading to more intense particle segregation. It is also observed that with 
increase in particle size or decrease in heat transfer coefficient, there is a reduction in 
particle rich region - this has aheady been discussed earlier in sections 4.1.2 and 4.1.6. It 
is noted that, the heat transfer coefficient and particle size are two important controlling 
parameters in deciding the particle segregation pattern. With increase in heat transfer 
coefficient, the heat removal rate increases leading to decrease in solidification time 
resulting in the increase in particle rich region. The final thickness of the particle rich 
region for AI-AI 2 O 3 varies from 4.524x10'^ to 2.44x10'^ m and 4.702x10'^ to 2.55x10'^ 
m for the values of heat transfer coefficient 2000 and 5000 W/m^/K respectively, as the 
particle size increases from 1 to 3 /jm. 

As solidification proceeds, thickness of the air-gap formed at the metal/mold 
interface increases due to the contraction of the casting as well as thermal expansion of 
the mold. This results in the reduction of the heat transfer coefficient at the interface with 
time. Thus the rate of solidification decreases as the solidification front progresses from 
the outer periphery of the casting towards its center - this is evident in Figs. 4.18 to 4.23. 
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Figure 4.18; Variation of volume fraction with time for 1 /m particle size 
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Figure 4.19; Variation of volume fraction with time for 2 jjm particle size 
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Figure 4.20: Variation of volume fraction with time for 3 /jm particle size 



Figure 4.21: Variation of volume fraction with time for 1 jjm particle size 
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Figure 4.22: Variation of volume fraction with time for 2 fjm particle size 
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Figure 4.23: Variation of volume fraction with time for 3 /Jm particle size 
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4.2 CALCULATION OF SOLIDIFICATION TIME: 


For various sets of operating conditions, the time required for complete 
solidification of the casting has been calculated. The results of these simulations are 
presented in the form of total solidification time. 

4.2.1 Effect of Initial Pouring Temperature: 

Figure 4.24 shows the effect of imtial pouring temperature on the complete 
solidification time for three different particle sizes, namely, Ipn, IjLon and Zfjm. The 
initial pouring temperatures taken in the simulations are 650, 730 and 800 °C. All other 
parameters are kept the same. It is observed that the total solidification time increases 
with increase in the initial pouring temperature. This is due to the fact that with increase 
in the initial pouring temperature the total heat content in the melt increases and thus 
more time is required to withdraw this extra heat. It is also seen that the total 
solidification time increases with increasing particle size. It has already been discussed in 
the earlier section that particle segregation becomes more intense with increase m the 
particle size. Since AI 2 O 3 particles are heat-resistant, increase in the volume firaction of 
these particles near the outer periphery of the casting decreases the thermal dififusivity in 
that region substantially. This results in decrease in heat removal rate, which, in turn, 
leads to increase in the solidification time. 

4.2.2 Effect of Initial Mold Temperature: 

Figure 4.25 is a plot between the initial mold temperature and solidification time 
for three different particle sizes, namely, l//m, Ifjm and 3jjm. The initial mold 
temperatures used in these simulations are 200, 250 and 300 °C. All other parameters 
remain the same. From the figure it is evident that the time required for complete 
solidification of the casting increases with increase m the initial mold temperature for all 
the three particle sizes. As the mitial mold temperature increases, the temperature 
gradient between the casting and the mold decreases reducing the heat transfer rate. This 
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4.2.3 Effect of heat transfer coefficient at the metal/mold interface: 

Figure 4.26 shows the effect of heat transfer coefficient at the metaiymold 
interface on the complete solidification time for three different particle sizes. The hi 
values used here are 1000, 2000 and 5000 W/m^/K. All other parameters are the same. It 
is observed that with increase in the heat transfer coefficient the complete solidification 
time decreases as the rate of heat removal becomes faster. 

4.2.4 Effect of Particulate Volume Fraction: 

Figure 4.27 shows the variation of total solidification time as a function of initial 
particulate volume fraction. As the initial volume firaction of the particulate increases, the 
effective thermal diffusivity of the composite decreases resulting in longer solidification 
time. 

4.3 TEMPERATURE DISTRIBUTION IN THE CASTING AND MOLD 
REGIONS: 

The temperature profiles in the casting as well as mold regions have been 
calculated by solving the governing differential equations with appropriate boundary 
conditions along with the energy balance equation discretized in Finite Volume 
Technique, which is described in the chapter 3. Figure 4.28 shows the variation of 
temperature as a function of time for both the casting and mold regions. It is seen that in 
the final temperature distribution, there is a temperature drop at the casting/graphite as 
well as at the graphite/steel mold interfaces. 

Figure 4.29 shows the effect of variation of heat transfer coefficient on 
temperature distribution in both the casting and mold regions. Here, the heat transfer 
coefficients used are 1000, 2000 and 5000 W/m'/K. All other parameters are fixed as 
earlier. It is observed that the temperatures at the inner and outer surfaces of the graphite 
mold increase with the increase in the heat transfer coefficient. This attributes to the fact 
that the higher heat transfer coefficient signifies better thermal contact between the 
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Figure 4.26: Effect of heat transfer co-efficient on solidification time for different particle 
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Figure 4.27: Effect of initial particiilate volume fraction on solidification time 
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Figure 4.29: Effect of heat transfer co-efficient on the temperature distribution for both 
the casting and mold regions 
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casting and the mold wall resulting in more intense heat transfer from the casting to the 
mold. On the other hand, the reduction in the heat transfer coefficient at the interface 
results in lesser heat transfer and therefore lower temperature at the interfaces. 

4.4 DISCUSSIONS: 

Siva Raju and Mehrotra [9] have presented a mathematical formulation based 
on one-dimensional heat transfer analysis incorporating variation in thermo-physical 
properties due to particle movement in the matrix. The present investigation is an attempt 
to further improve that analysis. Both Kang and Rohatgi [7] and Siva Raju and Mehrotra 
[9] in their models neglected the repulsive force term in their equations for force balance 
on particles. In the present investigations, the repulsive force, which a particle 
experiences when it approaches a solid wall/solidification front is taken into 
consideration and characterized by an appropriate expression. Again, in the work of Siva 
Raju and Mehrotra, variation of volume fraction of particles with time is considered by 
ensuring conservation of total mass of particles at all times. In their work, the volume of 
particulates in between two consecutive nodal points always remains same. So, m that 
case the length of the segment varies with time and so are the nodes as they depend on 
constant volume fraction of particulates. But in the present formulation the nodes are 
fixed for all time and the variation of volume fraction of particles in any particular 
segment with time is calculated by taking into account the particle movement with time. 
One of the important features of the present model is the use of finite volume technique as 
the solution procedure unlike the finite difference scheme used in the previous work [10]. 
The advantage of finite volume method over finite difference one is the automatic 
conservation of all physical quantities unlike the approximation of so in the later case. 

It is observed in sections 4.1.1 to 4.2.4 that the particle segregation depends on the 
rotational speed of the mold, size of the reinforcing material, relative density difference 
between the matrix and reinforcement material, initial pouring temperature of the liquid 
melt, pre-heatmg temperature of the mold, heat transfer coefficient at the casting/mold 
interface and the initial volume fraction of the particulates in the melt. The solidification 
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time gets influenced mainly by the heat transfer co-efficient at the casting/mold interface 
and initial temperatures of the liquid melt and mold. 

Formulation of the model is based on assumptions, several of which are only 
partially valid. For instance, the assumption of interface between the solid and liquid 
region to be planar is only approximate. Strictly speaking, the solidification front can be 
planar, cellular or dendritic depending upon the rate of solidification. It is very difficult to 
predict either cellular or dendritic morphology of the solidification front, hence the 
assumption of planar solidification front is made. 

The assumption of negligible thermal resistance between the particles and melt 
can also be questioned. In actual practice, there may be some solubility of solid 
particulates in the liquid matrix, or some chemical reaction between the two. In both 
cases, there may be a heat of reaction associated with the phenomena taking place at the 
surface of the particles. This heat may appear as the heat source at the particle/melt 
interface. Further, a chemical reaction at the surface may result in the formation of a 
product layer which offers thermal resistance to the heat flow, resulting in distortion of 
the temperature field aroimd the particle. This may, in any case, happen also due to the 
difference in thermal conductivities between the particle and the melt. 

The assumption of temperature invariant thermal properties is just for 
convenience only. If reliable expressions for thermal properties as a fimction of 
temperature are available, these can be easily incorporated in the model equations. 

Dependence of heat transfer coefficient at the casting/mold interface on the 
rotational speed is assumed to be negligible here. But it has been estabhshed that for 
aluminium, the heat transfer coefficient at the metal/mold interface increases with 
increase in the rotational speed [34]. The reason is that at a higher rotational speed, liquid 
metal puts a larger pressure on the solidified layer which, in turn, results in a better 
physical contact between the solidified layer and mold wall, improving the heat transfer 
coefficient at this interface. So, the effective heat transfer coefficient at the metal/mold 
interface can be found out using a suitable relation and this relation along with Eq. (2.14) 
will give a better estimate of variation of the effective heat transfer coefficient with time. 

For simplicity, initial position of the particles are considered at the nodal points. 
But it is assumed that particles are uniformly distributed in the melt. Uniformity of 
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Figure 4.30: Particle segregation pattern w.r.t. the number of grids 



particulate distribution can be described more effectively with larger nodal points. Higher 
number of nodes imply lesser grid size. But there is a limit beyond which the grid size 
can not be lowered since it has to be greater than the particle diameter. Figure 4.30 shows 
the variation of particle segregation pattern in the solidified melt with number of grids for 
AI-AI 2 O 3 system. It is observed in the figure that, with higher number of nodes, clustering 
of the particles is seen closer towards the outer surface of the casting and for 1000, 1200 
and 1400 number of grids, the segregation pattern assumes a similar shape and almost 
near to one another. But it is noted that, when the number of grids are 10, the volume 
firacticn of particles in the solidified composite becomes equal to that of the imtial melt. 
This is because, with a highly reduced number of grids, the distance between the particles 
increases. As a result of which particles have to travel farther distances to enter into the 
adjacent volume segment, which require a longer time. Here it can be inferred that, 
particles do not get proper time to move to the adjacent volume segment to form cluster 
and get trapped by the sohdification front when they are in the same segment as those of 
the initial ones. 

While frying to converge the time step in the program, several improvements have 
to be incorporated to achieve a much reduced convergence criteria ( lO"* ), to be ensured 
of the fact that the numerical technique is close to the exact solution. The improvements 
made were as follows. 

• To calculate the values of p, c and k for any interface of the control volume, 
harmonic mean of the adjacent nodal point values have been taken instead of 
arithmetic mean. 

• Throughout the program, all the calculations have been carried out with double 
precision to minimize the round-off error. 

• V alue of TT used here is 22/7 instead of 3 . 1 4. 
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CHAPTERS 


SUMMARY AND CONCLUDING REMARKS: 

The present investigation is carried out widi the objective of improving upon the 
already formulated mathematical model to describe the centrifiigal casting of metal 
matrix composites. The model aims to describe the temperature distributions throughout 
the casting and the mold regions and particle segregation pattern in the casting region. In 
the present investigation Finite Volume Technique has been employed in place of earlier 
used finite difference based solution scheme. In the force balance expression, one 
significant improvement in the model formulation is repulsive force term is incorporated 
for the particles that are at the vicinity of the solid/liquid interface. This modified 
expression is used to calculate the particle segregation pattern in the casting region. The 
influence of different process parameters such as rotational speed of the mold, particle 
size, relative density difference between the particle and the melt, initial mold 
temperature, pouring temperature, heat transfer coefficient at the casting/mold interface 
on the solidification time of the casting and on the segregation pattern of the particles is 
studied. 

The main conclusions of this investigations can be summarized as follows: 

• Thickness of the particle rich region in the composite decreases with increase in 
rotational speed, particle size, relative density difference between the particle and 
the melt, initial pouring temperature and initial mold temperature. 

• With decrease in the heat transfer coefficient at the casting/mold interface, the 
rate of heat transfer decreases leading to increase in the solidification time which, 
in turn, results in more intense segregation of solid particulates. 

• With increase in the initial volume fraction of solid particulates, both the 
solidification time as well as the final thickness of the particulate rich region 
increase. 

• Rate of solidification is very fast at the beginning, which is more than that of the 
particle velocity due to centrifugal force. Due to this, the particles adjacent to the 
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casting/mold interface do not get a chance to move towards the casting/mold 
interface in the cases of AI-AI2O3 and Al-SiC systems to form a cluster. As a result 
of which, the volume fraction of the particulates in the solidified composite is the 
same as that of the initial volume fraction in the melt which can be visualized by 
the horizontal portions of the figures in the cases otAl-AhOs and Al-SiC systems. 
So, in case of AI-AI2O3 and Al-SiC, the maximum clustering of particles is seen 
after a certain distance from the casting/mold interface. But in case of Al-Gr 
system, particles move towards the iimer periphery and form maximum 
segregation there. It is also seen that for 1 fan particle size in case of Al-Gr 
system, the volume fraction of the particles in the solidified composite, in the 
volume segment, just adjacent to that of the casting/mold interface remains same 
as that of the initial volume fraction in the melt unlDce to that of 2 and 3 pm 
particle size. 
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